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We consider the small separation asymptotic expansions of the Casimir interaction energy and 
the Casimir interaction force between two parallel cylinders. The leading order terms and the 
next-to-leading order terms are computed analytically. Four combinations of boundary conditions 
are considered, which are Dirichlet-Dirichlet (DD), Neumann-Neumann (NN), Dirichlet-Neumann 
(DN) and Neumann-Dirichlet (ND). For the case where one cylinder is inside another cylinder, the 
computations are shown in detail. In this case, we restrict our attention to the situation where the 
cylinders are strictly eccentric and the distance between the cylinders d is much smaller than the 
^— ' distance between the centers of the cylinders. The computations for the case where the two cylinders 

are exterior to each other can be done in the same way and we only present the results, which turn 
^ • up to be similar to the results for the case where one cylinder is inside another except for some 

^ ' changes of signs. In all the scenarios we consider, the leading order terms are of order d~'^ ^"^ and 

. they agree completely with the proximity force approximations. The results for the next-to-leading 

order terms are new. In the limiting case where the radius of the larger cylinder approaches infinity, 
the well-known results for the cylinder-plate configuration with DD or NN boundary conditions are 
recovered. 



I. INTRODUCTION 



Recently, there has been an increasing interest in the Casimir effect from both the theoretical and the experimental 
sides [l| . Before the turn of this century, the theoretical studies on the exact Casimir effect were mostly restricted to 
ps) . simple geometries such as parallel plates, spherical shells and cylindrical shells. However, these geometric configura- 
J> ' tions pose a certain degree of difficulty for the experimental verification of the Casimir effect, such as the difficulty 
in achieving parallelism. As a result, experimentalists favor other configurations, especially the sphere-plane config- 
\l uration. However, in the last century, one has to rely on the proximity force approximation to estimate the Casimir 
^ ^ \ force between such configurations, which hampers the determination of the experimental accuracies. To circumvent 
. • the problem, intensive activities have been carried out to research for theoretical methods to determine the Casimir 
25 interactions between two or several objects beyond the accuracies afforded by the proximity force approximations. 
In the last decade, a number of methods have been developed, which include the semi-classical approximation 0,11], 
I ' the optical pat h method the worldline approach 0"1O|, the functional determinant or the multiple scattering 

• • , method pl]42Cll |. and the exact mode summation method 21,[22|. Using the multiple scattering method, one can in 



>^ ' principle write down a functional for the Casimir interaction energy between two or several objects. Nevertheless, it 
Jv>( , should be mentioned that the mode summation approach can also lead to the same result, although this latter method 

; I ■ has only been applied to the configuration of two eccentric cylinders. 

^ In principle, using the various methods mentioned above, one can compute the magnitude of the Casimir force 
numerically. However, the accuracy is always subjected to the computing capacity of the computer, especially when 
the separation between the objects is small, which is of more interest for comparison to experiments. On the other 
hand, to determine the dependence of the Casimir force on various parameters of a configuration, one usually has to 
assume a priori the form of the dependence on the parameters and determine the coefficients that fit best for that 
particular form. These coefficients are subjected to numerical errors and it is not easy to justify the accuracy of the 
results. Therefore there is a call for analytically computing the asymptotic expansion of the Casimir force when the 
separation between the objects is small. For the cylinder-plate and the sphere-plate configurations, the first corrections 
to the proximity force approximations have been computed analytically in pTl . [23l - [25| . Although the configurations 
of two cylinders and two spheres are among the popular ones whose exact Casimir interaction energies have been 
derived using the multiple scattering approach or the mode summation method [13, [3412) l26l - [28| , analytical studies 
on the corrections to the proximity force approximations for these configurations are still lacking. The purpose of this 
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article is to address this problem for the configuration of two cylinders. The results would be compared to the results 
for th e cy linder-plate configuration which is a special case where the radius of one of the cylinders approaches infinity. 



In [28|, |29| , some experimental setups have been proposed to measure the Casimir interaction force between two 
eccentric cylinders. It was argued that the cylindrical configuration has some of the experimental advantages of 
both the the parallel planes and the sphere-plane configurations. It can lead to favorable conditions to search for 
extra-gravitational forces in the micrometer range and for the observation of finite-temperature corrections. These 
latter subjects have been explored for the cylinder-plane configuration in j30l - [33| . In view of this, the cylindrical 
configuration is becoming increasingly important in the study of the Casimir effect. It is timely to do an analytical 
study on the strength of the Casimir force when the separation between the cylinders is small. 

In this article, we consider the case where one cylinder is inside another, and the case where two cylinders are 
outside each other. The functional determinant representations of the exact Casimir interaction energies have been 
derived in several places, such as [20l - [23 , l33 - [36j for the case of one cylinder inside another, and [l|, HO, HE Hll for 
the case of two cylinders exterior to each other. Using this formulas, we compute analytically the first correction to 
the proximity force approximations. In the case of one cylinder inside another, the roles of the two cylinders are not 
symmetrical. We explain the computations in detail for this case. The computations for the case of two cylinders 
outside each other can be done in the same way and we only present the results. 
Throughout the paper, we use the units with h = c = 1. 

II. PROXIMITY FORCE APPROXIMATION OF THE CASIMIR INTERACTION BETWEEN TWO 

PARALLEL ECCENTRIC CYLINDERS 



As shown in Fig. [1] we consider two parallel cylinders of length L and radii a and b respectively. The cylinder of 
radius a lies inside the cylinder of radius b. Denote by S the separation between the centers of the cylinders, and d the 
distance between the cylinders. Clearly, we have S = b — a — d. In this article, we assume exclusively S > 0, i.e., the 
cyfinders are not concentric. The concentric case has been considered in [H, [13, 111] ■ For the boundary conditions on 
the cylinders, we impose Dirichlet or Neumann boundary conditions. In the following, we will denote the boundary 
conditions by XY, where X = D (Dirichlet) or N (Neumann) is the boundary condition on the cylinder of radius a, 
and Y = D or N for the cylinder of radius b. 

In this section, we use the proximity force approximation (PFA) to obtain the leading term of the Casimir interaction 
force between the two cylinders when d ^ (5 <C 6 — a. Notice that the radii a and b are fixed parameters, whereas the 
distance between the cylinders d is a variable. 

First, recall that the Casimir interaction force per unit area on two parallel plates separated by a distance H is 
given by 





FIG. 1: The cross section of two eccentric cylinders. 
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if the two plates are both imposed with Dirichlet boundary conditions or Neumann boundary conditions. If one 
plate is imposed with Dirichlet boundary condition and one plate is imposed with Neumann boundary condition, the 
Casimir interaction force density is —7/8 times of ([T]), i.e., 

II _ 7 TT^ _ 77r2 

cas( ) - 8 480iJ4 - 3840i?4 ' ^ ' 

Using PFA, we have to integrate the Casimir energy density J^Q^^iH) over the area of one of the cylinders, with H 
being the distance from a point of the integrated cylinder to the other cylinder. Here we choose to integrate over the 
cylinder of radius b. The integration over the length of the cylinder is trivial. Using polar coordinates, we find that 
the shortest distance from the point with parameter 6 (see Fig. 1) to the cylinder of radius a is 



y/{bcos0-6)^ + (6 sin 61)2 - a = y/b^ + 6^ - 26(5 cos 6* - a. 
Then for the DD or the NN case, the PFA for the Casimir interaction force between the cylinders is 



FEl"---^! ^d9. 



240 



(V62 + ^2_ 26(5 cos e - a) 



Making a change of variables 



we find that 

62 + (52 - {du + a)2 



\/62 + (52 - 26(5cos6i 
d 



cos 6 —- 



2b6 



JAb^6^ - [62 + (52 - (du + a)2]2 ^d(u - l)(2a + du + d) [(26 - a - d)^ - (du + a)2 
sin^?=-^^ — 



26(5 26(5 



d[du + a) 



and thus 



2(b-a)-d 

^PFA TT bL I << d(du + a) 2 ^ 

Cas -~^^J^ ^i^i —^^^ _ -^)(2a + du + d) [(26 - a - d)2 - {du Ta)^ 



In the limit d 0, we find that the leading order term of the PFA is 



^PFA TT^V^L 1 

240^/2(6 -a)d5 A u'^y/ll^ 
TT^VabL 



(3) 



768v/2(6 - a)(i2 



For two perfectly conducting eccentric cylinders, the proximity force approximation to the Casimir interaction force 
is twice that of ([3]), corresponding to the sum of the TE (which is equal to NN) and the TM (which is equal to DD) 
contributions. This has been obtained in 29]. 

For the DN or the ND case, one just has to multiply ([3]) by —7/8 which gives 



^PFA Tn'^VabL 



" ^ . (4) 

6144^2(6 -a)di 



III. THE FORMULA FOR THE EXACT CASIMIR INTERACTION ENERGY BETWEEN THE 

PARALLEL ECCENTRIC CYLINDERS 



The exact Casimir interaction energy between two eccentric cylinders has been derived in [211, l22|, |3J| using mode 
summation approach and in [20,, ,35i . ,36i] using scattering or functional determinant method. The Casimir interaction 
energy can be represented as 

Ecs^j^ ^Trln(l-M(0)d^, (5) 
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where M is an cx) x oo matrix with elements Mm„, — oo < m, n < oo 



P— — 00 



^--^^^ KM) m) '^-"-^'^^'^-^'^^^ 

oo , (6) 

In these formulas, Ii,(z) and are the modified Bessel functions of first and second kinds respectively. Recall that 

the dependence on the distance d between the cylinders is encoded in the variable S — b — a — d. From ^ , we notice 
that the terms dependent on d, Ip-m{SS.)Ip-niSS,), is identical for the four combinations of boundary conditions. The 
factors Ip-miSS,) and /p_„(i5^) actually come from translation formulas. 

Expanding the logarithm and taking the trace in the formula ([S]), we find that the Casimir interaction energy can 
be evaluated as 

s— ^ Jo — — oo js— — oo 

Alternatively, one can also use the identity Tr In = In det to write the Casimir interaction energy ^ as 



T POO 

Ecs^— ^lndet(l-M(0)de 
47r Jo 



(8) 



In principle, the magnitude of the Casimir interaction energy can be evaluated with the help of a computer. When the 
separation between the cylinders d is much larger than the radius a of the smaller cylinder, one can in fact determine 
analytically the dominating term of the Casimir interaction energy from ([7]), since in this case, the dominating term 
is the term with s = and jo = 0. In the opposite limit where d <C a, 6, which is of more experimental interest, 
even the numerical computation of the magnitude of the Casimir interaction energy or the Casimir interaction force 
poses a great challenge since one has to take the matrix M with larger size in ([8|) for convergence. Some numerical 
results have been discussed in [s^ for the special cases of quasiconcentric cylinders (where S <^ a), concentric cylinders 
(where ^ = 0) and the limiting case of a cylinder in front of a plane (where b ~> oo). In the following, we will use 
analytical approach to compute the leading order term and the first order correction term of the Casimir interaction 
energy and the Casimir interaction force when d <C a, 6. In addition, we also assume that d <C 6 — a. However, we do 
not make any assumption about the relative sizes of S and a. As mentioned above, we assume exclusively that S ^ 0. 
In fact, d <C 6 — o implies d <C (5 and S ^ b — a. The method we use is similar to that used in [111, li^l for a cylinder 
in front of a plane and in [2^ [23, [4l[ for a sphere in front of a plane. Compare to the case of a cylinder in front of a 
plane discussed in [ill . |40| , the main complication here is the summation over p that appears in the matrix elements 

IV. ASYMPTOTIC BEHAVIOR OF THE CASIMIR INTERACTION ENERGY AND THE CASIMIR 
INTERACTION FORCE OF THE ECCENTRIC CYLINDERS AT SMALL SEPARATION 

Define the dimensionless parameter 

a 

" = I 

— a 



and the dimensionless variable 



d 



Let 

/3 = a + 1 = 



b ~ a 

b 



b — a 
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We are interested in the asymptotic behaviors of the Casimir interaction energy and the Casimir interaction force 
when £ <C 1. Making a change of variables f = ^/{b — a) in (O, we find that the Casimir interaction energy can be 
rewritten as 



Eca 



/ c. — r\ ' "^0 „■ — o — -n^ — _, , — _rv-^ 



47r(6- „ , 

S — ^ jo — — CXD js— — OOpo — — OC Ps— — 00 



where A„„;p(w) = Bmn;pi^)Tmn;p{uj) , 

ODD / N _ /n(acj) gp(/3a^) rNN / N /^(a^) ^^l/?^) 

DN . . _ /n(«^) ND . . _ (10) 

T'm„;j,(w) ==/p_„((l - e)a;)/p_„((l - e)uj). 

To find the first two leading terms of the Casimir interaction energy, one can replace the summations by integrations, 
i.e., 

OO „oo pOO pOO s s s 

(2s+2) times 

with the understanding that j'^+i = jo- 

Using Debye asymptotic expansions of modified Bessel functions one finds that A^-^^^.^^^luj) has an expansion 
of the form 



A 



(w) = ± — exp(^z^i7](a;i) + i'2V{^2) - 2(^377(^3) + v^-qiuji) + v^ri{u)^)j 



I'l 1^2 It's 

Z, 



where 



1^1 1/2 1^3 1^4 / ' 



1^4 = - i^l, = - 1^2, 

acj aw /3a; (1 — e)w (1 — e)a; 

wi = — , W2 = — , i:j3 = — , 1^4 = , = ; 

1^1 V2 v\ 



(12) 



r]{z) ^ Vl + z2 + In i(z) 



iDD _ 4DN _ I ( 1+W;^ 



J^2J^4J^5 V(l +'-'2)(l +'-'4)(l +^^5)/ ' 



iNN _ /iND I V2 ( 1 + ^^2 



J/1 1^41^5 V(l +^^l)(l +'^4)(1 +wi)/ ' 

2;c(0 — ""1(0 (o^ •2c(i) = ^i(O) if the cylinder of radius c is imposed with Dirichlet boundary condition (or Neumann 
boundary condition), with 

-it , ^ It^ - 9t 

Ul (t) = , V-] (t) = . 

The plus or minus sign in (jl2l) depends on the boundary conditions. For DD or NN boundary conditions, we have 
the plus sign; whereas for DN or ND boundary conditions, we have the minus sign. 
When the variables z/i , 1/2, 1^3, w, e vary, the function 

vir]{uJi) + V2V{^2) - 2iy3i]{uJ3) + v^rfluji) + vsvi^s) 
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is always nonpositive. It achieves the maximum value of when vi = 1/2, e = and 

a a 

In that case, 

This suggests that in ([9]), we can rename jo as m, and introduce new variables ni, . . . , n<,, qo, . . . ,qs so that 

ji = m + fii, 1 < « < s, 
/3 , 

Pi = — (2to + n,j + nj+i) + g.,, < i < s. 
la 

By a further substitution 



we find that 



^^"^ 3=0""+^ Jo Jo J-00 J-oo,,=o i^o ^=l ^ 



(13) 



(2s+l) times 



with the understanding that no — Ug+i — 0. 
With 



/3 , 

z/i=TO + ni, n2 = m + ni+i, V3 = —{2ni + ni + n,+i) + q,, 

Za 

and treating ni,ni^i,qi and e as perturbed variables, one can easily deduce that viri{uii) + ^211(^2) — '^ty^rj^ujs) + 
Vi-q^LJi) + v^rj^uj^) has an perturbative expansion of the form 

C30 C30 

i^iri{uJi) + V2ri{L02) - 2z^377(w3) + Vi-q{uji) + v^riiio^) X! 51 £^'m^^''Gkj{ni,ni+i,qi)'Hkj{T), (14) 

k=Q j=Q 

where Gkj{ni,ni^i,qi) is a homogeneous polynomial of degree k in ni,ni^i,qi. The terms with (A:,j) = (0,0) and 
= (1,0) are identically zero. The leading terms are —971;, where 

SH, = +l^(n,-n,+i)V^V- 

ar Am mp 



(fcj) = (0,l) (fcj) = (2,0) 

comes from the terms with {k,j) — (0, 1) and (fc, j) ~ (2, 0). This implies that the leading contribution to the Casimir 
interaction energy comes from m ^ e^^ , rii, ri^+i, qi ^ ^fm ^ . Hence, the (fc, j) term in (ITi)) is of order g-J+'^/^-i. 
To keep everything up to order e, we need the terms in ([14]) with j + k/2 < 2. With the help of a computer, we find 
that up to terms of order e, 

where 21^ and 03; are terms of order ^/e and e given respectively by 

^^=^^\-^^ + + — (-.-n,+i) + g j-er[2q, + 

^ T^{iT^ -I) ( a'^{a'^ + 2,a + i)qt [a + 2){n, + n,+^)q^ a^qj {{a^ + a){ni - m+if + [n, + m+if) 

JJi — ^ — ' 



to3 \ 12;33 6/3 

a^qt{ni + ni+i)(ni - Ui+iY P 



192 

er(l-T2)/ 2 / ^ (a^(nj -nj+i)2 + (nj +nj+i)2)\ ^2^^ 

m \ ia I a 
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In the same way, for A, one finds that 



A (1 + e, + 

m 



where <Li and S)^ are respectively terms of order y/e and e given by 

TO 

OTO^ V ^ ^ / 

in the DD case. The DN case is the same as the DD case, and the NN or ND case can be obtained from (|15p by 
interchanging rii and rii+i. Finally, one can check that the leading term of 

Ui 1^2 V5 

is of order e. Denote by the leading term. In the DD case, it is given by 

DD ^ (l + a + a^)T(5r2~3) 
12m/3 



whereas for the NN, DN and ND case, we have 



to/3 ' 
ar(T2 - 1) 
/3m ' 

* * to 

Notice that does not depend on Ui^rii^i and Collecting the terms, we find that 

A™+„^,™+„.+,;p^(w(t)) - ± ^- exp(-m, + 21, + (1 + + D,) (l + S^,) 

ZTTTO V / 

where 25^ and Sji are respectively terms of order ^/e and e given by 

2 

Substituting into the Casimir interaction energy (|13p . we find that up to first order correction term. 



Cas 



°^ ( 1 \v(s+l) pi poo pOO pOO ( s \ 



(2s4-l) times 

s s—1 s s \ s s 

1 + ^ 0i + ^ ^ ©^©j + XI n n dnidnidr. 

i=0 i=0 j=i+l 4=0 / 1=0 j=l 

Here x = for DD or NN boundary conditions, and x — ^ for DN or ND boundary conditions. The leading order 
term of the Casimir interaction energy is 



po _ 



00 / iw(s+l) rl poo poo poo I ^ \ ^ * 

(2s+l) times 
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The integrations of the term 

s 



4=0 



of order ^ over ni^rii^i^qi give zero since this term is odd in either qi^Ui or n^+i. Thus the next-to-leading order 
term is 



Cas 2^+2 



5+1 r °° ( I'lxCs+i) /■! r°° r°° f " \ 



(2s+l) times ^J^y^j 

' s — 1 s 

^ ^ ©,;0j + XI -^^^ I n '^^ n driidmdT, 

^i=0 j=i+l i=0 j i=0 i=l 

which is of order e smaller than the leading order term. 

Let us first consider the leading order term (|16p . It is straightforward to compute. Integrating over qi,0 < i < s, 
first, we have 

oo poo ( s ^IJl^ \ s i±i ^ r>£+i 



exp 



mB 
1=0 ^ 



times 

For the integrations over rii, 1 < i < s, since [ll|: 

— ? 

3 



^K-n,+i)2 = 2('ni-in2) + | fr^2 - ^^s) + . . . + ^^^n^, (18) 



j=0 

and 



y g "^o^^dx = y'^ for any a;o, 



(19) 



integrating in the order ni — > ri2 — > . . . — )• ris gives 



\ »=0 / i=l 



j32T2^S + 1 



Therefore, 



. - (_l)x(s+i) .1 3 / 2(s-M)em\ , 

-SSas - - 2^ \ ,,3 / ^ M m2 exp ^— ^ dmdT 

47r2a2 t:^ s + 15 Jo Jo V / 



Saiy^L ^ (-l)x(^+i) /-i 



64\/27ra2£5 ^ (s + 1) 



(20) 



Using the fact that 

1 _ TT^ (~1)^ _ 7 TT^ _ TtT^ 

^(s + l)4~90' ^ (s + 1)4 ~ 8 90 ~ 720' 

s=0 ^ ' s=0 ^ ' 

we find that for DD or NN boundary conditions, the leading order term of the Casimir interaction energy is 



1920v/2(6-a)di' 
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whereas for DN or ND boundary conditions, 



E'^c^^ = / 22 

lbim^2{b - a)di 



For the Casimir interaction force, one then obtains 



for DD or NN boundary conditions, and 



FLs "'^"^ , (23) 

^ 768^2(6 -a)d5 



^^^as = '"'^"^ , (24) 
^ 6144^/2(6 - a)di 



for DN or ND boundary conditions. The leading terms of the Casimir interaction forces and agree completely 
with that derive from proximity force approximations. 

For the next-to- leading order term (jl7p . we integrate first with respect to qi using the formulas 



oo 



0, if i is odd, 

q'e-^'i'dq^l r(2+i) (25) 



2 

A 2 



if j is even. 



A straightforward computation gives 



2 



s times 
' s— 1 s 



i=0 j=i+l i=0 / 1=1 



(26) 



where 



®^ / exp 2^ 

The term i^i can be decomposed into two parts: 

where S'i is as before and is independent of ni and n^+i. The explicit expressions for and ^ are given in Appendix 
El They are functions of and n^+i. Let 



s times 



s times 



and 



i=0 1 = j=i+l 
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so that 



^^Cas = - 2^ \/, ,,3 ' i exp ( ^— ^ ) S^dmdr. (27) 



47r2a2 ^ (s + 1)^ 



The integrals Zf and J^j are Gaussian and can be integrated in exactly the same way as in [ll|, which we explain 
in Appendix [Bl For the DD case, we find that 



3a2/3 

{[-7{s + 1)2 + 3a + 2]t2 + 4(s + 1)2 + ^2 - a - l) 



16/3m(s + l) 

For the DN case, since = ©f^ and 

pm 

we find that 

^s,DN ^ gs,DD _ a(s + l)r(T2 - 1) 

For the NN case or the ND case, ©f ^ = ©f ° and J?f ^ = j^nd ^^^^ obtained from 6°° and J?P° respectively by 
interchanging and As explained in our previous work [40], these imply that J^j^^ — Jij^^ — Js'-^/^s-i ^'^^ 

xs,NN _xs,DD ^NN t;DD _ rs,DD ''"(''"^ ~ 1) 
-'i --'s-z - -^s-i ^ ^^^^ ' 

js.ND ^js,DD ^ g^ND _ ^DD _ js.DD ^ ''"('''^ ~ 1) 
1 s 1 1 1 s 1 7"7^ 

Hence, 



g..NN^^s.DD^ (i + lMl^J)^ (28) 

to/3 

^.,ND ^^s,DD ^ + IM^' - 1) ^ (29) 

m 

We notice that the change in due to the change of boundary conditions only comes from the term involving ui (t) 
and vi{t) in the Debye asymptotic expansions of the modified Bessel functions and their derivatives. 
Now, it is easy to perform the integration over m and r in ([27]). We find that in the DD case, 

1 DD _ _ 3v^L ^ 1 ( (24a+16)r^ + 7a^-a-3 \ 



3VabL ^ 1 7 7 \ 



64V27r(6-a)idi ^ (s + l)Hl2 36a/3 
TT^V^L / 7 7(6 -a) 



1920^2(6 -a)di V 12(6 -a) 36a6 
In the NN case, dH]) implies that 



3\/^L ^ 1 



VPL 1 /-i^^i r^l^^ ( 2(s + 1)£to \ (s + l)T(r2-l) 



64V^7r(6-a)idi ^ + 1)^ 

TT^V^L / 7 



1920^2(6- a)d5 V 12(6 -a 



7 


7 


_^( 


12^ 


36a/3 


" 7 


40 ■ 


6-a\ 


36 


" 3^ 


a6 y 



(31) 
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Here we have used the fact that 



TT 



For DN and ND boundary conditions, the summation over s is alternating in sign. Using 



we have 



E, 



1,DN 
Cas 



3VabL 



i-iy 



TT 

12' 



■E 



64V2^(& - a)idi ^ (s + 1)^ V 12 36a 



7 7 

+ 



7 



160 1 



15360v/2(6 -a)(ii V 12(6 -a) 
Finally, for the ND case, it is easy to find that 



7 b-a 



36 ah 2l7r2 6 



(32) 



1,ND 

Cas 



7Tr^V^L 



15360i/2(5-a)rfi V 12(6 -a) 



7 b-a 
36 a6 



160 1 

2l7r2 a 



(33) 



Combining with the leading terms, we find that the asymptotic expansions of the Casimir interaction energies are 
given, up to the first corrections, by 



^Cas 



TT^VabL 



1920^2(6 - a)S 



abL 



TP 

E 



DN 

Cas 



1920^2(6 -a)rfi 
7TT^yML 



1+d 

1 + d 



ciND 
■'^Cas 



15360^2(6- a)di 
7t:^ VabL 



1 



15360^/2(6- a)d2 
For the Casimir interaction forces, we then obtain 

'!r^\/abL 



-^Cas 



-f^Cas 



Cas 



rjiND 

Cas 



768^2(6 -a)cii 

■jT^VabL 
768^2(6 - a)di 

7n^VabL 



1 



6144a/2(6 - a)d5 
77r3V^L 



6144v'2(6 - a)<i5 



1 + d 

l + d 

l + d 



12(6 -o) 
7 

12(6 -a) 
7 

12(6 -a) ^ 
7 

12(6 -a) "* 



7 

20(6 -a) 
7 



+ 



7(6 - a) 



36a6 



+ 



40^ 

36 37r2 



6- a 



ab 

160 r 



7(6 -a) 

36a6 ^ 2l7r2 h 

7(6 -a) 160 1 



(34) 



20(6 -a) 
7 



+ 



36a6 



7(6 - g) 
60a6 

7 i 
60 ~ 71 



2l7r2a 



6 — a 
ab 



+ 



20(6 - a) 60a6 

7 7(6 - a) 



7(6 -a) 32 1 

^7^6 



(35) 



20(6 -a) 



60a6 



32 1 
77r2 a 



From these, it is easy to see that if the smaller cylinder is imposed with Dirichlet boundary conditions (the DD or DN 
case), then the proximity force approximation underestimates the strength of the Casimir interaction force. However, 
if the smaller cylinder is imposed with Neumann boundary conditions, then the proximity force approximation may 
overestimates or underestimates the strength of the Casimir interaction force depending on the relative sizes of the 
two cylinders. Numerically, 



20(6 - a) 
7 



+ 



7 8 
60 ~ 7r2 

7(6 - a) 



20(6 - a) 60a6 



b-a __ 0.6939(6 - 0.4985a) (6 - 2.0059a) 
a6 a6(6 — a) ' 

32 1 _ 0.3465(6+ 0.1815a)(6- 1.8549a) 
777^ a ab{b — a) 
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Therefore there is a critical ratio of h/ a over which the proximity force approximation overestimates the strength of 
the Casimir interaction force. 

An interesting limiting case to study is the cylinder-plate configuration which can be achieved by taking the limit 
6 — )■ oo. In this case, we find from psp that the first two leading terms of the Casimir interaction force are given by 



F 



CP,DD 



Cas 



F. 



CP,NN 



Cas 



F 



CP,DN 



Cas 



F 



CP,ND 



Cas 



768\/2di 
768V2di 
' QlAAy/2di 
'6144^/2(^i 



1 + — - 



7_d 
60a 
7_ 
60 
J_d 
60a ^ 
7 

60 " 



32 
7^ 



The results for the DD case and the NN case have been obtained in [ll| 
results for the DN and the ND case have not been obtained before. 



(36) 



However, to the best of our knowledge, the 



V. TWO PARALLEL CYLINDERS EXTERIOR TO EACH OTHER 



In this section, we consider the case where two cylinders are parallel and exterior to each other (see Fig. [2]). The 
results for this case can be obtained in the same ways as the case of one cylinder inside another cylinder which we 
consider above. Now, the distance between the centers of the cylinders S is related to the radii a and b of the cylinders 
and the distance between the cylinders dhy6 = a + b — d. The proximity force approximation shows that at small 




FIG. 2: The cross section of two cylinders exterior to each other. 

separation (i.e., d <^ a, fo), the leading term of the Casimir interacting force between the cylinders is given by 

- 'V^bL 

^ 768^2{a + b)di 

for DD and NN boundary conditions. For DN and ND boundary conditions, we have 



^PFA ^_!^^^^^^_ ^38^ 
6144^2(0 + 5)^5' ^ ' 

The exact Casimir interaction energy between two cylinders which are exterior to each other has been derived using 
scattering approach or functional determinant method in P, [20l [ssl [36j . It can still be written in the form 



Ec.s^^ ^Trln(l-M(0)d^, 



(39) 
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where now the elements of the infinite matrix M are 



'KM) 



■ V — — OC ' 



(40) 



E 



'KM) ,^^K,m 



The computations of the leading order term and the next-to-leading order correction term of the Casimir interaction 
energy go in parallel to the former case. Let 

d 



a + b 

be the dimensionless variable. Then as e — >■ 0+, we find that 

TT^VabL 



ciDD 



77.NN 



■'^Cas 



1920^y2{a + b)di 

1920y/2{a + b)di 



ND 
Cas 



15360^/2(0 + 



l + d 
l + d 

+ d 



12(a-|-6) 
7 

12{a + b) 
7 



15360^2(a -|- b)di 
For the Casimir interaction force, we then obtain 

-K^VabL 



l + d 



12(a + b) 
7 

12{a + b) 



+ 



7{a + b)^ 
36ab 


-) 




(7 40 \ a + 


b 


V36 37r2 ) ah 


7(a + b) 


160 1' 




36a5 


2l7r2 b 


+ 


7{a + b) 
36a6 


160 1" 
2l7r2a 


+ 



^Cas 



^Cas 



Cas 



rjiND 

Cas 



768v'2(a -t- b)d2 



768 V2(a + b)di 

7Tr^VabL 



l + d 
l + d 



7 



6144 v/2(a + b)d-2 
7'K^VabL 



6UAy/2{a + b)d2 



1 



l + d 



20{a+b) 
7 

~ 20{a + b) 
7 

'20ia+b) ^ 
7 



+ 



+ 



7{a + b) 



(41) 



60a6 

7_ _ 
60 



ab 



7{a + b) _ 32 1 
60a6 Ttt^I 

7{a + b) 



+ 



(42) 



20(a-F&) 60a6 



32 1 
77r2 a 



Notice that for the DD and the NN cases, there is a complete symmetry between the parameters a and b. On the 
other hand, the ND case can be obtained from the DN case by interchanging the parameters a and b. These are 
expected since in the present case, the two cyUnders arc on equal footing. 

Numerically, 7/60 — S/tt^ = —0.6939. Therefore, it is easy to conclude that in the NN case, the proximity force 
approximation overestimates the strength of the Casimir interaction force. For the DD boundary conditions. 



7 



20(a + b) mob 



7(a + b) 7(a^ -ab + b^) 

+ . = — TT^ > 0. 



60a6(a -|- b) 



Therefore, the proximity force approximation underestimates the Casimir interaction force. For the ND boundary 
conditions. 



7 



7{a + b) 32 1 
Wa + bj ^ 60ab Tir^a 



0.3465(6-1- 1.8549a) (6- 0.1815a) 

ab{a + b) 
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Therefore, we find that the proximity force approximation may underestimate or overestimate the strength of the 
Casimir interaction force depending on the ratio of the two radii of the cyhnders. 

It is interesting to compare the results for two cylinders exterior to each other (j4T|l and (|42p with the results for one 
cylinder inside another cylinder and pSI) . Notice that they have similar coefficients up to the changes of signs. 
This can be considered as an analogy of the result obtained in [l^, where the exact closed forms for the Casimir 
energy of two weakly coupled dielectric cylinders were derived and it was shown that the analytic results for one 
cylinder inside the other can be obtained as an analytic continuation of the result for two cylinders exterior to each 
other. 

Taking the limit b — ^ oo, we again recover the configuration of a cylinder in front of a plate. It is easy to check that 
taking the 6 ^ oo limits of dH]) give ([M)) . 



VI. CONCLUSION 



In this article, we have computed analytically the asymptotic expansion of the Casimir interaction force between 
two cylinders, with one inside the another, or both outside each other. We compute the leading order term and 
the next-to-leading order term. Different combinations of Dirichlet (D) and Neumann (N) boundary conditions are 
discussed. The results read as 



^Cas 



-^Cas 



r^DN 

Cas 



abL 



768y/2{b Ta)di 



±- 



7 



abL 



768y/2{bTa)S 
7n^VabL 



± 



F, 



ND 

Cas 



6144^2(6 Ta)rf^ 
7'K^VabL 



1 



6144^/2(6 Ta)rf^ 



± 



± 



"20(6Ta) 
7 

20(6 T a) 
7 

20(6 T a) ^ 
7 



60a6 



6=F a 



60 7r2 ) ah 

7(6 T a) ^ 32 l' 
60a6 7^6 

7(6 T a) 32 1 

20(6 =F a) ^ 60a6 7^ a 



(43) 



For the terms ± or =p, the sign on the top is for the case where one cylinder is inside another, and the sign at the 
bottom is for the case where the two cylinders are exterior to each other. It is observed that for each case, the leading 
term of the Casimir interaction force agrees with that derived using the proximity force approximation. The proximity 
force approximation may underestimate or overestimate the magnitude of the Casimir interaction force depending on 
the boundary conditions and the ratio of the radii. The special b ^ oo limiting case which gives the results for the 
cylinder-plate configuration is discussed. It is found that the 6 — > oo limits of the asymptotic expansions for the DD 
and NN cases reproduce the well-known results for the corresponding asymptotic expansions for the cylinder-plane 
configuration. 

Although wc only consider the Dirichlet and Neumann boundary conditions in this article, it is easy to obtain from 
(|43p the asymptotic expansions for perfectly conducting or infinitely permeable cylinders. More specifically, if both 
the cylinders are perfectly conducting or both are infinitely permeable, one takes the sum of the results for the DD 
and the NN cases. If one cylinder is perfectly conducting and one is infinitely permeable, then one takes the sum of 
the results for the DN and the ND cases. 



Appendix A: Explicit expressions for &i and 

For &, 



■'1 1 



gDD _ gDN _ r^(nt - 3ni+i) ^ /?r'^(nj + n.i+i){ni - Uj+iY _ ^ r(n^ -|- n^+i) 
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and and 25^^^ are obtained from this by interchanging and rii+i. For 

= j^^^^^ Yq^2 ((a +a)(n, -n,+i) +(n, +n,+i) ) 

r3(3r2 - 1)13 



192m3 



4^ - + + ) + ^(1 - - ) + ^ 



^2 



8rn2 

5(a + 2)V3 a(a + 2)T4 _ 2 £(a + 2)T2 Sr^K + n.+i)^ /?r^(n| - n^+i)^ 

48m/? 16m2 2a 32m2 32m3 

+ - ""'^'^ - - + 128;;^3 



Bm^a 2a'^ Am Srn^ 



(rii-ni+i) + 



— — -\- £ — . 

Am? 8m^ ma a 

Kf^ and are obtained from this by interchanging rii and ni+i. 



Appendix B: The integrals /| and Jfj 

For /q, using the identity 

X^/ n2 o / 1 V 2 \^ s+1 , ,^ , 

2_^[ni~ni+i) =2lns--ns_il + - I n^-i - -ns_2 I +-.-H n^, (Bl) 



we can integrate in the order n2 ni. The integrations over rt^, ... ,712 only require the formula (|19p . 

whereas the integration over ni needs the formulas (j25p . 

For /| , we use (flS)) and integrate in the order rii —s- . 

For /f , 1 < i < s — 1, using the identity 

^(rii - rij+i)^ =2 ^ni - ^"2^ + . . . + ^-y^ ^Ji^ - — ^"-i+i^ + 2 ^ri^ - ^"s-i^ + • ■ • 

s — i — 1 \^ s + 1 



'^i+a : —ni+i ) + ^. , ^ 



s— i— 1\ s — i / (i + l){s — i) 

we can first integrate in the order ni — > . . . — n,:_i —>■... —> and then make a change of variables 



X ^ n.i - ———ni+i, y^m+i. 
I + 1 



The integrations over x and y can be performed using 

It is interesting to remark that although the cases of Iq and /| have to be considered separately, it turns out that 
by formally substituting i = and i — s into the formula obtained for where 1 < i < s — 1, the results agree with 
Jg and /| respectively. 

For Jq^ , using the identity (|B1[) , we can first integrate in the order —>...—> 77,3 . Then make a change of variables 

s-1 

x = ni, y = n2 ni, 



and integrate over x and y using ([25 
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For Jqjj, using the identity 

ini-n^+i) =2lm2--m3l +...H — ^ I ms_2 — 1 H ^ I m-s-i — j-("s-ni) 



1\^ s-2/ s-3 V s-l f s-2 



1 



2 V S-3V s-2 y S-2V s-1 

1 V s+l 2 
f^s H 



where rrii ~ rii — ni, 2 < i < s — 1, we can first integrate in the order TO2 — > . . . — > ?7t.s-i- Then make a change of 
variables 

1 

x = ni Hs, y^Us, 

s 

and integrate over x and y. 

For Jqj , where 2 < < s — 1, using the identity 

2 



i=0 



J - 2 V ' 3 

( 1 \^ s — 2 ( s — j — 1 

+ 2i^ns- -n,_i j + . . . + J— —J - _ ^ n^+i 

if 1 ^ \ ^ + 1 2 , - i + 1 -"^ - .? 

ni + — ; —n^ H — n^+i 



where = — ni, 2 < i < — 1, we can first integrate in the order 7712 — > — > . . . — > nj+2- Then 
make a change of variables 

1 s — 

X = ni rrij, y = rij, z = nj+i -nj, 

J s — J + 1 

and integrate over x,y and z. 

For Js_x s, using the identity (|18p. we can first integrate in the order 7ii — > . . . ns-2- Then make a change of 
variables 

s-1 

X = Us^i y 

s 

and integrate over x and y. 

For Jfs, 1 < i < s — 2, using the identity 



/ 1 
^(rii - n.i+if =2 ( ni - -712 



1=0 



2 . / . -, X 2 

I — 1 



z — 1 \ i 

( 1 \^ s — z — 1/ s — i — 2 

2 TOs-i - 7:'7is-2 + • • • + ^ ^ : ririi+i - n^) 

\ 2 / s — I — 2 \ s — I — 1 

/ \ 2 / 

i + 1 f i \ s + 1 9 s — i I 1 



l + l J [l + l)[s—l) S— 2— 1 * 

where rrij = rij — n^, i + 2<j<s — l,we can first integrate in the order rzi —> — mg-i —>■...—> m,j4.2. 

Then make a change of variables 

i 1 

n-i+i, y = n,+i, z = Ug r^i+i, 



and integrate over x, y and z. 

For J.fi^i, 1 < « < s — 2, using the identity 



j=0 



1 \ 2 •11/ • N 2 

1 \ 2 + 1 / I 



^(rii - Ui+i f =2 i ni - -n2 j + . . . H — ( - — — n^+i 



i \ i + l 

1 \^ s — i f s — i — 1 \^ s + 1 2 
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we can first integrate in the order m ^ . . . — )• n,_i ni+3. Then make a change of variables 

X = Ui - ——-rit+i, y^rit+i, z = ni+2 —rii+i 

I + I s — I 

and integrate over x, y and z. 

Finally, for J^J, where l<i<s~3,i + 2<j<s — 1, using the identity 

1 \ 2 • , 1 / • \ 2 



j=0 



y^(rti - rij+i)^ =2 ( m - -n2 ) + • • • H — ( - ——ni+i 



+ 2 ris - -ns_i + . . . H — rij+i r— r%- 

V 2 ; s-j \ s-j+1 

( 1 \^ 3 — ( J — « — 2 

+ 2 ( mi+2 - 2'^i+3 ) + ■ • • + j_i_2 y^^-^ ~ ~ 



where nik = Uk — fii+i, « + 2 < fc < j — 1, we can first integrate in the order ni ^ . . . — >■ nj_i ^ — ^ 

mj+2 —>...—>■ rrij-i. Then make a change of variables 

i z ~\~ 1 s — J 

x = ni- —-rrii+i, y = n^+i —nj, z = nj, w = nj+i _ nj, 

and integrate over x, y, z and w. 

As in the case of If, by formally substituting (i, j) = (0, 1), (0, s), (0, j), (s — 1, s), (i, s), (i, i + 1) into the expressions 
obtained for J?- where l<i<s — 3,i + 2<j<s — 1, one obtains respectively -'oi' -^Os' •^Oj' -^s-i.s' -^is' "^ti+i- 
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